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\ We update the constraints on new-physics contributions to AF = 2 processes from 

' the generalized unitarity triangle analysis, including the most recent experimental 

■ 

■ developments. Based on these constraints, we derive upper bounds on the coefficients 



of the most general AF = 2 effective Hamiltonian. These upper bounds can be 



H ' translated into lower bounds on the scale of new physics that contributes to these 

low-energy effective interactions. We point out that, due to the enhancement in the 
renormalization group evolution and in the matrix elements, the coefficients of non- 
standard operators are much more constrained than the coefficient of the operator 
present in the Standard Model. Therefore, the scale of new physics in models that 
generate new AF = 2 operators, such as next-to-minimal flavour violation, has to 
be much higher than the scale of minimal flavour violation, and it most probably 
lies beyond the reach of direct searches at the LHC. 



I. INTRODUCTION 



Starting from the pioneering measurements of the K'^ — mass difference Am^ and of 
the CP- violating parameter ek, continuing with the precision measurements of the Bd — Bd 
mixing parameters Attlb^ and sin 2/5 and with the recent determination of the Bg — Bg 
oscillation frequency Am^^ and the first bounds on the mixing phase — 2/?^, until the very 
recent evidence of — Z)° mixing, AF = 2 processes have always provided some of the 
most stringent constraints on New Physics (NP). 

For example, it has been known for more than a quarter of century that supersymmetric 
extensions of the Standard Model (SM) with generic flavour structures are strongly con- 
strained by — mixing and CP violation jl|. The constraints from — mixing 
are particularly stringent for models that generate transitions between quarks of different 
chiralities 0, 3, 4]. More recently, it has been shown that another source of enhancement 
of chirality-breaking transitions lies in the QCD corrections [s^, now known at the Next-to- 
Leading Order (NLO) [6 



9. 



Previous phenomenological analyses of AF = 2 processes in supersymmetry |8|, l9|] were 
affected by a large uncertainty due to the SM contribution, since no determination of the 
Cabibbo-Kobayashi-Maskawa lOj (CKM) CP-violating phase was available in the presence 
of NP. A breakthrough was possible with the advent of B factories and the measurement 
of time-dependent CP asymmetries in B decays, allowing for a simultaneous determination 
of the CKM parameters and of the NP contributions to AF = 2 processes in the and 
Q, Furthermore, the Tevatron experiments have provided the first 



Bd sectors 
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measurement of Am and the first bounds on the phase of Bs — Bs mixing. Combining all 
these ingredients, we can now determine allowed ranges for all NP AF = 2 amplitudes in 
the down-quark sector. 

To complete the picture, the recent evidence of — mixing allows to constrain NP 
contributions to the AC = 2 amplitude HQ. 

Our aim in this work is to consider the most general effective Hamiltonian for AF = 2 
processes {Ti.^^"^) and to translate the experimental constraints into allowed ranges for the 
Wilson coefficients of Ti.^^^. These coefficients in general have the form 

a(A) = 4^ (1) 



where Fj is a function of the (complex) NP flavour couplings, Li is a loop factor that is present 



in models with no tree-level Flavour Changing Neutral Currents (FCNC) , and A is the scale 
of NP, i. e. the typical mass of the new particles mediating AF = 2 transitions. For a generic 
strongly-interacting theory with arbitrary flavour structure, one expects ~ Lj ~ 1 so that 
the allowed range for each of the Ci{A) can be immediately translated into a lower bound on 



A. Specific assumptions on the flavour structure of NP, for example Minimal 



Next-t o- Minimal 



3, 
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18| or 



19j Flavour Violation (MFV or NMFV), correspond to particular choices 



of the Fi functions, as detailed below. 



Our study is analogous to the operator analysis of electroweak precision observables 20| , 
but it provides much more stringent bounds on models with non-minimal flavour violation. 
In particular, we find that the scale of heavy particles mediating tree-level FCNC in models 
of NMFV must lie above ~ 60 TeV, making them undetectable at the LHC. This bound 
applies for instance to the Kaluza-Klein excitations of gauge bosons in a large class of models 
with (warped) extra dimensions 2l| . Flavour physics remains the main avenue to probe such 
extensions of the SM. 

The paper is organised as follows. In Sec. [TTlwe briefly discuss the experimental novelties 
considered in our analysis. In Sec. IIIII we present updated results for the analysis of the 
Unitarity Triangle (UT) in the presence of NP, including the model-independent constraints 



on AF = 2 processes, following closely our previous analyses 



11 



I'M . In Sec. [IV]we discuss 



the structure of Ti.^^"^, the definition of the models we consider and the method used to 
constrain the Wilson coefficients. In Sec. |V]we present our results for the Wilson coefficients 
and for the scale of NP. Conclusions are drawn in Sec. I VI I 



II. EXPERIMENTAL INPUT 



We use the same experimental input as Ref. [12[, updated after the Winter '07 conferences. 
We collect all the numbers used throughout this pager in Tables U and HIl We include the 
following novelties: the most recent result for Am, 
decays [^23,] and the dimuon charge asymmetry from D0 



22l |. the semileptonic asymmetry in B 



24 and CDF [251], the 



measurement of the Bg lifetime from flavour-specific final states 26|, the determination of 
AFs/Ps from the time-integrated angular analysis of Bs J/'4"P decays by CDF the 
three-dimensional constraint on F^, AF^, and the phase 0s of the Bs-Bg mixing amplitude 
from the time-dependent angular analysis of Bg —>■ J/ip4> decays by D0 [28^. 



Parameter 




Value 


Gaussian (a) Uniform 

(half-width) 


A 




0.2258 


0.0014 


iFefel X 103(excl.) 




39.1 


0.6 1.7 


iKfel X 103(incl.) 




41.7 


0.7 


iKfel X 104(excl.) 




34 


4 


iKfel X 104(mcl.) 




43.1 


3.9 


Arrid (ps~^) 




0.507 


0.005 






17.77 


0.12 


ei^ X 10^ 




2.232 


0.007 


sin 2/3 


see Winter 


'07 analysis at http://www.utfit.org 


cos 2/3 


see 


Winter 


'07 analysis at http://www.utfit.org 


a 


see 


Winter 


'07 analysis at http://www.utfit.org 


7 


see 


Winter 


'07 analysis at http://www.utfit.org 


2/3 + 7 


see 


Winter 


'07 analysis at http : // www . utf it . org 


(GeV) 




161.2 


1.7 


asiMz) 




0.119 


0.003 


TBs (ps) 




1.39 


0.12 


(ps) 




1.454 


0.040 


Ad 
^SL 




-0.0005 


0.0056 


4^ 
^SL 




0.0245 


0.0196 






-0.0043 


0.0030 






0.009 


0.037 


Ar,/r, 




0.65 


0.33 



(/)« [rad] -0.79 0.56 (Ps-tBs corr. 0.727 

TB, (ps) 1.49 0.08 TB^-ATs/Ts corr. -0.172 

ATs/Ts 0.17 0.09 ATs/Ts-cps corr. -0.188 



TABLE I: Values of the experimental input used in our analysis. The Gaussian and the flat 
contributions to the uncertainty are given in the third and fourth columns respectively (for details 
on the statistical treatment see 29|]). See text for details. 



Parameter Value Gaussian (cj) Uniform 

(half-width) 



Fd (MeV) 201 
FB^VB'siMeY) 262 



e = 

Bk 

rrib (GeV) 
me (GeV) 
Ri 

R2 

R4 

R5 



BP 
B? 

bF 

BF 
BF 



B 



Bi 



bF 

Bi 
Bf 
Bi 



1.23 
0.79 
4.21 
1.3 
1 

-12.9 
3.98 
20.8 
5.2 
0.865 
0.82 
1.07 
1.08 
1.455 
0.88 
0.82 
1.02 
1.15 
1.99 



3 
35 
0.06 
0.04 
0.08 
0.1 

3.0 
0.89 
4.4 
1.2 
0.02 
0.03 
0.05 
0.02 
0.03 
0.04 
0.03 
0.06 
0.03 
0.04 



17 



0.08 



0.015 
0.01 
0.08 
0.02 
0.075 
0.10 
0.09 
0.13 
0.13 
0.24 



TABLE II: Values of the hadronic parameters used in our analysis. The Gaussian and the flat 
contributions to the uncertainty are given in the third and fourth columns respectively (for details 



on the statistical treatment see 



29I]). See the text for details. 



The use of AFo/r.,, from the time-integrated angular analysis of Bg J /'4"P decays, is 
described in Ref. 12] ■ In this paper, we only use the CDF measurement as input, since the 
D0 analysis is now superseded by the new time-dependent study. The latter provides the 



first direct constraint on the Bg-Bg mixing phase, but also a simultaneous bound on AF^ 
and r^. We implemented the full 3x3 correlation matrix. The time- dependent analysis 
determines the Bg-Bg mixing phase with a four-fold ambiguity. ^ First of all, the Bg mesons 
are untagged, so the analysis is not directly sensitive to sin0s, resulting in the ambiguity 
(0s, cos 2) ^ (— 0s, — cos(5i,2), where 61^2 represent the strong phase differences between 
the transverse polarization and the other ones. Second, at fixed sign of cos5i^2; there is the 
ambiguity {(f)g, ATg) ^ {(f)g + tt, — AF^). Concerning the strong phases 6i, there is a two-fold 
ambiguity corresponding to 6i 71 — 6i. The two experimental determinations are roughly 



(^2 ~ 0, (5i ~ TT and 61 ~ 0, 62 



TC. 



corresponds to the first choice [30|, 



Si 
B 



0, (^2 ~ 



33 



34 



31 



n the literature it is often found that factorization 



321]. However, we find that factorization predicts 



35| . This result is also compatible with the BaBar measurement in 
which can be related to Bg —>■ J/'^cf) using 5'f/(3) and neglecting singlet 



contributions.^ However, waiting for future, more sophisticated experimental analyses which 
could resolve this ambiguity, we prefer to be conservative and keep the four-fold ambiguity 
in our analysis. 



The use of Ams was already discussed in Ref. 12|. The only difference with respect 
to that is the update of the experimental inputs: we now use the improved measurement 
by CDF 2^, and we take tb^ only from the study of Bg decays to CP eigenstates 37|. 



The value of tb^ obtained from Bg decaying to flavour-specific final states, using a single 



exponential in the fit, is related to the values of Fs and AF^ by the relation 



r. 



FS 



1 1 + 



Bs 



F 



^ 1 



V 2r» 



2 ' 



(2) 



which provides an independent constraint on AFs/F^. We compute AFs and using 
eq. (7) of Ref. (recalling that = 2(1 — \q/p\g). Following Ref. 39], we use the value 
of recently presented by D0 2J] and CDF 25|, in the form 



A 1^1^ 



fdXdO + fsXsO 



(3) 



1 Notice that the definition used by D0 is the one of Ref. [301, namely = 2/5^ = 2 a.rg{- {VtsV^f,) / (VcsV;^))) 

in the SM. Notice also that in the arXiv version of Ref. 28] the definition of 0s is unclear. 
^ In the first version of this manuscri pt, we stated that factorization disagreed with SU{3), based on the 
factorization prediction in Refs. 



30 



31 



with fd = 0.397 ± 0.010, fg = 0.107 ± 0.011, Xqo = iXg + Xg)/2- Xq and Xg are computed 
using equations (3)- (5) of Ref. [3].^ 

Finally, concerning — mixing, we use as input the results for the NP amplitude 
obtained in Ref. jl^ combining the experimental information from Refs. 40|. 



III. UT ANALYSIS AND CONSTRAINTS ON NP 



The contribution of NP to AF = 2 transitions can be parameterized in a model- 



independent way as the ratio of the full (SM+NP) amp^ 
we can define the parameters Cn and 0d (g = d, s) as 



itude to the SM one. In this way. 



4l|: 



{Bq 


rrfuU 


Bq) 




Bq) 



and write all the measured observables as a function of these parameters and the SM ones 
(p, f/, and addition al p arameters such as masses, form factors, and decay constants). Details 



are given in Refs. H 



121]. In a similar way, one can write 



Im[(iro|ifSM|^o)] ' ^^^K Re[(KO|^SM|^o)] • 

Concerning Amx, to be conservative, we add to the short-distance contribution a possible 
long-distance one that varies with a uniform distribution between zero and the experimental 
value of ArriK- 

We perform a global analysis using the method of Ref. I29I and determine simultaneously 
p, f], Cb^, (pBg, CArriK usiug flat a-priori distributions for these parameters. The 

resulting probability density function (p.d.f.) in the p — fj plane is shown in Fig. [H Only 
a small region close to the result of the SM fit survives. The mirror solution in the third 
quadrant is suppressed down to about 5% probability by the measurements of and A'^j^. 
The results for p and fj reported in Tab. Illll are at a level of accuracy comparable to the SM 



fit 42[, so that the SM contribution to FCNC processes in the presence of arbitrary NP is 
bound to lie very close to the results of the SM in the absence of NP. This result represents 
a major improvement in the study of FCNC processes beyond the SM, and opens up the 
possibility of precision studies of flavour processes in the presence of NP. 



To combine the CDF and D0 measurements, we have converted the value for A defined in Ref. [241 into 
a value for Ag[^. 



FIG. 1: Determination of p and 77 from the NP generalized analysis. 68% and 95% probability 
regions for p and f\ are shown, together with the 2 a contours given by the tree-level determination 
of I Kb I and 7. 



The constraining power of this analysis is evident in the results for the NP parameters 
given in Tab. IIIII and shown in Fig. [21 Cornpared to our previous analysis in Ref. 12| . 



and to similar analyses in the literature 



39, 



431], we see that the additional experimental 



input discussed above improves considerably the determination of the phase of the Bg — B 



mixing amplitude. The fourfold ambiguity inherent in the untagged analysis of Ref. [28] 
is somewhat reduced by the measurements of and ^sl, which prefer negative values 



of (pBs-'^ Ref- 22] recently claimed a 2a deviation from zero in 0^^, taking the sign of 



cos5i.2 from factorization. We confirm that, with the same assumptions of Ref. 32] on 
strong phases,^ the deviation from zero of (pBs slightly exceeds 2a. Without assuming strong 



^ With respect to Ref. (32| , we find that the inclusion of Ag^ has a weaker impact in reducing the ambiguity 
coming from several small differences in the analysis (theoretical assumptions on NP in Agj^, presence of 
NP in penguin contributions to ^4^'^, inclusion of the CDF measurement of ^gL, etc.). 

^ We find that factorization gives Si and 62 ^ tt, resolving the ambiguity of the DO untagged analysis 



Parameter 


68% Probability 


95% Probability 




1.05 ±0.34 


[0.53,2.05] 




-3.4 ±2.2 


[-8.3,1.6] 




1.11 ±0.32 


[0.63, 2.07] 




(-69 ± 14) U (-20 ± 14) U (20 ± 5) U (72 ± 8) 


[-86,-46] U [-43,35] U [56,87] 


CAmfc 


0.93 ±0.32 


[0.51,2.07] 


CeK 


0.92 ±0.14 


[0.66,1.31] 


P 


0.140 ±0.046 


[0.049, 0.258] 


V 


0.384 ± 0.035 


[0.304, 0.460] 


an 


86 ± 7 


[73, 104] 




24.3 ± 2.0 


[19.9,28.6] 




70 ± 7 


[53, 83] 


ReAt[10-^] 


-31.3 ± 1.8 


[-35, -27] 


ImAt[10-5] 


14.8 ± 1.3 


[12.3,17.3] 


\Vub\[W-''] 


3.91 ±0.28 


[3.36, 4.44] 


\Vcb\[W-^] 


4.09 ± 0.05 


[3.98,4.18] 




8.7 ±0.5 


[7.6,9.6] 


\Vtd/Vts\ 


0.217 ±0.012 


[0.187,0.240] 


Rb 


0.413 ± 0.030 


[0.354,0.471] 


Rt 


0.943 ± 0.048 


[0.819,1.037] 


sin 2/3 


0.748 ± 0.044 


[0.643,0.841] 


sin 2Ps 


0.0409 ± 0.0038 


[0.0321,0.0494] 



TABLE III: Determination of NP and UT parameters from the UT fit. 
phases from factorization, we find a deviation of (pBs from zero of 1.4 a, of the same size of 



the deviation found experimentally by Ref. 28]. 



It is important to stress that the information contained in these constraints does not rely 



in favour of (ps ~ 0.79 for positive AF, while Ref. 84] uses i5i ~ tt, ^2 ~ and (f>s ~ —0.79. However, this 
sign difference in (j)s is compensated by the fact that (j)s as defined in Ref. [34] should be compared to — 
as measured by DO. 




FIG. 2: Constraints on (f)B^ vs. Cs^, (pBs vs. Cb,, and Cf^j^ vs C^mK from the NP generalized 
analysis. 



on any specific model for NP. The list of applications in NP phenomenology is rich. For 
instance, restricting to the case of SUSY models, the two-dimensional constraint of the b s 
{b d) sector can be translated into a limit on the mass-insertion complex parameters 
(5^3), using the NLO supersymmetric expression of the Bq-Bn mixing amplitude [3]. This 
bound, combined with the constraint from b —>■ s decays 4J], allows to obtain the best 
available information on the off-diagonal terms of the squark mass matrix [45]. However, 
in the following, rather than considering an explicit model for NP, we perform a general 
analysis based on the most general Hamiltonian for AF = 2 processes. 



IV. NP CONTRIBUTIONS TO AF = 2 PROCESSES 



The most general effective Hamiltonians for AF = 2 processes beyond the SM have the 
following form: 

ntn'='=j:QQf+j:QQf (6) 
1=1 1=1 

5 3 _ _ 



i=l 1=1 



1=1 i=l 

where q = d{s) for Bd(s) — Bd[s) mixing and 

r)<li<lj _ ^ 

_ p,a -/3 



Qr^ = <ifi^^LQ?i^?L , (7) 



/n-Ji* _ ^ ri"' ri^ 

V4 — IjRliL'djL^iR 

f^lili _ ri^ p^P rP" 

^5 — HjRQiLHjLHiR 



Here gR^i, = Pr^lQ, with Pr^^ = (1 ± 75)/2, and a and /3 are colour indices. The operators 
2^3 are obtained from the *3 by the exchange L ^ R. In the following we only discuss 
the operators Qi as the results for Qi^2,3 apply to (5i,2,3 as well. 

The NLO anomalous dimension matrix has been computed in [6|. We use the 
Regularisation-Independent anomalous dimension matrix in the Landau gauge (also known 
as RI-MOM), since this scheme is used in lattice QCD calculations of the matrix elements 
with non-perturbative renormalization. 

The Ci{A) are obtained by integrating out all new particles simultaneously at the NP scale 
A.^ We then have to evolve the coefficients down to the hadronic scales fib = = 4.6 GeV 
{rrih = rrihifi = rrif,) is the RI-MOM mass) for bottom mesons, fic = 2.8 GeV for charmed 
mesons, and fix = 2 GeV for Kaons, which are the renormalisation scales of the operators 



used in lattice computations for the matrix elements 46l. |47|. 



^ Clearly, without knowing the masses of new particles, one cannot fix the scale A of the matching. However, 
an iterative procedure quickly converges thanks to the very slow running of as at high scales. 



We give here an analytic formula for the contribution to the Bg — Bg mixing amplitudes 
induced by a given NP scale coefficient Ci(A), denoted by {Bg\Ti.c{i\Bg)i, as a function of 
«s(A): 

{B,\n^/='\Bgh = E E (6f ^ + r/cf v""^ Q{A) {BglQ^W , (8) 



j=l r=l 



where rj = as{A)/as{mt), the magic numbers aj, bj~'^^ and c^^'^' and the matrix elements can 



be found in eqs. (10) and (12) of Ref. respectively. The values of the Bj^ parameters can 
be found in Table [Tll A similar formula holds for — Z)° mixing, with the parameters -Bf 
given in Table [III and the following magic numbers: 

tti = (0.286,-0.692,0.787,-1.143,0.143) 



(11)^ 

i 


(0.837,0,0,0,0), 


(11) 

Ci = 


(-0.016,0,0,0,0), 


(22) _ 
i 


(0,2.163,0.012,0,0), 


„{22) _ 


(0,-0.20,-0.002,0,0), 


(23) _ 
i 


(0,-0.567,0.176,0,0), 


„(23) _ 


(0,-0.016,0.006,0,0), 


(32) _ 
i 


(0,-0.032,0.031,0,0), 


^(32) _ 


(0,0.004, -0.010,0,0), 


(33) _ 
i 


(0,0.008,0.474,0,0), 




(0,0.000,0.025,0,0), 


(44) _ 
i 


(0,0,0,3.63,0), 


^(44) _ 


(0,0,0,-0.56,0.006), 


(45) _ 
i 


(0,0,0,1.21,-0.19), 


„{45) _ 


(0,0,0,-0.29,-0.006), 


(54) _ 
i 


(0,0,0,0.14,0), 


„(54) _ 


(0,0,0,-0.019,-0.016) 


(55) _ 
i 


(0,0,0,0.045,0.839), 


J55) _ 


(0,0,0,-0.009,0.018). 



bf = (0, -0.032, 0.031, 0, 0), cP = (0, 0.004, -0.010, 0, 0), (9) 

All other magic numbers vanish. Finally, for — mixing we obtain 

mnti='\K% = E E (^f ^ rf^ cm (irO|gf I^O) , (10) 

jr = l r=l 

where now the magic numbers can be found in eq. (2.7) of Ref. j^. We use the values in 
Table HI] for the ratios Rr of the matrix elements of the NP operators Q^*^ over the SM one. 
These values correspond to the average of the results in Ref. (46^ applying a scaling factor 
to the errors to take into account the spread of the available results. 

To obtain the p.d.f. for the Wilson coefficients at the NP scale A, we switch on one 
coefficient at a time in each sector and calculate its value from the result of the NP analysis 
presented in sec. IIIII 

As we discussed in eq. ([T]), the connection between the Cj(A) and the NP scale A depends 
on the general properties of the NP model, and in particular on the flavour structure of the 



Assuming strongly interacting new particles, we have from eq. ([T]) with Lj = 1 

A = \/|. (11) 
us now discuss four notable examples: 

• In the case of MFV with one Higgs doublet or two Higgs doublets with small or 
moderate tan/5, we have Fi = Fsm and Fj^i = 0, where Fsm is the combination of 
CKM matrix elements appearing in the top-quark mediated SM mixing amplitude, 
namely {VtgV^l)'^ for Bg — Bg mixing and {VidV^*)"^ for Ek- ^itlk and — mixing 
do not give significant constraints in this scenario due to the presence of long-distance 
contributions. 

• In the case of MFV at large tan /?, we have this additional contribution to Bg — Bg 

C,(A)= (°° + °-);f° + °'' A,A,FsM, (12) 



mixing 



where Xb^g represent the corresponding Yukawa couplings, ao,i,2 are tan /3-enhanced 
loop factors of 0{1) and A represents the NP scale corresponding to the non-standard 
Higgs bosons. 



In the case of NMFV, we have \Fi\ = Fsm with an arbitrary phase [19|] (following 
Ref. 181], for Aitlk and — mixing we take Fsm = l^d^sP)- This condition is 
realized in models in which right-handed currents also contribute to FCNC processes, 
but with the same hierarchical structure in the mixing angles as in the SM left-handed 
currents. Given the order-of-magnitude equalities m^/mb ~ \Vtd\, m^/mfe ~ \Vts\, 
bounds obtained in this scenario are also of interest for extra-dimensional models with 
FCNC couplings suppressed linearly with quark masses [2l|. Clearly, given the QCD 



and, for — K'^ mixing, chiral enhancement of NP operators, the constraints on the 
NP scale are much stronger for NMFV than for MFV, as shown explicitly in the next 
section. 

For arbitrary NP flavour structures, we expect \Fi\ ~ 1 with arbitrary phase. In this 
case, the constraints on the NP scale are much tighter due to the absence of the CKM 
suppression in the NP contributions. 



V. RESULTS 



In this Section, we present the results obtained for the four scenarios described above. 
In deriving the lower bounds on the NP scale A, we assume Lj = 1, corresponding to 
strongly-interacting and/or tree-level NP. Two other interesting possibilities are given by 
loop- mediated NP contributions proportional to or a^. The first case corresponds for 
example to gluino exchange in the MSSM. The second case applies to all models of SM- 
like loop-mediated weak interactions. To obtain the lower bound on A for loop-mediated 
contributions, one simply multiplies the bounds we quote in the following by as{A) ~ 0.1 or 
by aw ~ 0.03. 



Let us first consider MFV models and update our results presented in Ref . [UJ, [12] . In 
practice, the most convenient strategy in this case is to fit the shift in the Inami-Lim top- 
quark function entering Bj^, and mixing. We fit for this shift using the experimental 
measurements of Am^, Am^ and ex-, after determining the parameters of the CKM matrix 



with the universal unitarity triangle analysis ITj.^ We obtain the following lower bounds at 
95% probability: 

A > 5.5 TeV (small tan (3) , (13) 
A > 5.1 TeV (large tan/3) . (14) 

The bound for large tan/3 comes from contributions proportional to the same operator 
present in the SM. 

As mentioned above, at very large tan/? additional contributions to Ci{h) can be gener- 
ated by Higgs exchange. From these contributions, we obtain the following lower bound on 
the scale A, which in this case is the mass of non-standard Higgs bosons: 

Mh > 5 ^{ao + ai)(ao + a^) TeV . (15) 

In any given model, one can specify the value of the Oj couplings and of tan/5 to obtain a 
lower bound on the non-standard Higgs mass. The bound we obtained is in agreement with 



Ref. 



181 ]. taking into account the present experimental information. If a non-standard Higgs 



^ With respect to the original proposal of Ref. [l7|, we do not use the ratio Arris / Arrid in the fit in order 
to allow for Higgs-mediated contributions affecting Am^ at very large tan /3. 



boson is seen at hadron colliders, this implies an upper bound on the couplings and/or 
tan p. 

In Fig. [3] we present the allowed regions in the ReC*-ImC* planes for the sector, while 
in Figs, nun] we show the allowed regions in the AbsC*-ArgC* planes for the D^, Bd and 
Bs sectors. All coefficients are given in GeV~^. From these allowed regions we obtain the 
95% probability regions for reported in the second column of Tab. IIV[ This result is 
completely model-independent. 

Assuming strongly interacting and/or tree- level NP contributions with generic flavour 
structure {i.e. Li = \Fi\ = 1), we can translate the upper bounds on Ci into the lower 
bounds on the NP scale A reported in the third column of Tab. llVi As anticipated above, 
we see that in the K'^ sector all bounds from non-standard operators are one order of 
magnitude stronger than the bound from the SM operator, due to the chiral enhancement. 
In addition, operator Q4 has the strongest Renormalization Group (RG) enhancement. In 
the D^, Bd and Bg sectors, the chiral enhancement is absent, but the RG enhancement is still 
effective. The overall constraint on the NP scale A comes from ImCj^ and reads, for strongly 
interacting and/or tree-level NP, as loop mediated or aw loop mediated respectively: 

A,^™ > 2.4 ■ 10^ TeV, A^™ > 2.4 ■ 10^ TeV, A^™ > 8 ■ 10=^ TeV. (16) 

Assuming strongly interacting and/or tree- level NP contributions with NMFV flavour 
structure {i.e. Lj = 1 and = |-Fsm|), we can translate the upper bounds on Ci into 
the lower bounds on the NP scale A reported in the fourth column of Tab. IIVI The flavour 
structure of NMFV models implies that the bounds from the four sectors are all comparable, 
the strongest one being obtained from ImC^ (barring, as always, accidental cancellations): 
A^^v > Q2 TeV, A^f > g 3 TeV, A^^^^ > 3 TeV. (17) 

Let us now comment on the possibility of direct detection of NP at LHC, given the 
bounds we obtained. Clearly, a loop suppression is needed in all scenarios to obtain NP 
scales that can be reached at the LHC. For NMFV models, an aw loop suppression might 
not be sufficient, since the resulting NP scale is 2 TeV. Of course, if there is an accidental 
suppression of the NP contribution to ex, the scale for weak loop contributions might be as 
low as 0.5 TeV. The general model is out of reach even for aw (or stronger) loop suppression. 
For MFV models at large values of tan (3, stringent constraints on the mass of the non- 
standard Higgs bosons can be obtained. These particles may or may not be detectable at 




FIG. 3: Allowed ranges in the ReC|^-ImC]^ planes in GeV ^. Light (dark) regions correspond to 
95% (68%) probability regions. 






FIG. 4: Allowed ranges in the AbsCI^-ArgC}) planes in GeV ^. Light (dark) regions correspond 
to 95% (68%) probability regions. 





FIG. 5: Allowed ranges in the AbsC^^-ArgC^^ planes in GeV ^. Light (dark) regions correspond 
to 95% (68%) probability regions. 





FIG. 6: Allowed ranges in the AbsC^^-ArgC^^ planes in GeV ^. Light (dark) regions correspond 
to 95% (68%) probability regions. 



Parameter 


95% allowed range 


Lower limit on A (TeV) 


Lower limit on A (TeV) 




(GeV-2) 


for arbitrary NP 


for NMFV 


ReC^ 


[-9.6,9.6] • 10-13 


1.0 


103 


0.35 


ReC| 


[-1.8, 1.9] • 10-1^ 


7.3 


103 


2.0 


ReC| 


[-6.0,5.6] • 10-1^ 


4.1 


103 


1.1 


ReC^ 


[-3.6,3.6] • 10-15 


17- 


103 


4.0 


ReC|. 


[-1.0, 1.0] • 10-1^ 


10- 


103 


2.4 


ImC^ 


[-4.4,2.8] • 10-15 


1.5 


lO'' 


5.6 


ImCl 


[-5.1,9.3] • 10-1'^ 


10- 


10^ 


28 


ImC| 


[-3.1,1.7] • 10-16 


5.7 


10^ 


19 


ImC^ 


[-1.8,0.9] • lO-i'^ 


24- 


10^ 


62 


ImC|, 


[-5.2,2.8] • 10-1'^ 


14- 


10"^ 


37 




< 7.2 • 10-13 


1.2 


103 


0.40 




< 1.6 • 10-13 


2.5 


103 


0.82 




< 3.9 • 10-12 


0.51 


• 103 


0.17 




< 4.8 • 10-1"^ 


4.6 


103 


1.5 


1^^1,1 


< 4.8 • 10-13 


1.4 


103 


0.47 


nj 


< 2.3 • 10-11 


0.21 


•103 


1.8 


nj 


< 7.2 • 10-13 


1.2 


103 


8.5 


nj 


< 2.8 • 10-12 


0.60 


• 103 


4.4 


nj 


< 2.1 • 10-13 


2.2 


103 


14 


nj 


< 6.0 • 10-13 


1.3 


103 


8.8 


nj 


< 1.1 • 10-9 


30 


1.3 


nj 


< 5.6 • 10-11 


130 


4.6 




< 2.1 • 10-1° 


70 


2.4 




< 1.6 • 10-11 


250 


7.9 


\ck\ 


< 4.5 • 10-11 


150 


4.9 



TABLE IV: 95% probability range for C(A) and the corresponding lower bounds on the NP scale 
A for arbitrary NP flavour structure and for NMFV. See the text for details. 



the LHC depending on the actual value of tan (3. Finally, the reader should keep in mind the 
possibility of accidental cancellations among the contribution of different operators, which 
might weaken the bounds we obtained. 

VI. CONCLUSIONS 

We have presented bounds on the NP scale A obtained from an operator analysis of 
AF = 2 processes, using the most recent experimental measurements, the NLO formulae for 
the RG evolution and the Lattice QCD results for the matrix elements. We have considered 
four scenarios: MFV at small tan /?, MFV at large tan /?, NMFV and general NP with 
arbitrary flavour structure. The lower bounds on the scale A of strongly-interacting NP for 
NMFV and general NP scenarios (barring accidental cancellations) are reported in Fig. [71 
Taking the most stringent bound for each scenario, we obtain the bounds given in Table |Vl 




FIG. 7: Summary of the 95% probability lower bound on the NP scale A for strongly-interacting 
NP in NMFV (left) and general NP (right) scenarios. 

We conclude that any model with strongly interacting NP and/or tree- level contributions 
is beyond the reach of direct searches at the LHC. Flavour and CP violation remain the 
main tool to constrain (or detect) such NP models. Weakly-interacting extensions of the 
SM can be accessible at the LHC provided that they enjoy a MFV-like suppression of 
AF = 2 processes, or at least a NMFV-like suppression with an additional depletion of the 
NP contribution to ex- 



Scenario 


strong/tree 


Os loop 


aw loop 


MFV (small tan /3) 


5.5 


0.5 


0.2 


MFV (large tan f3) 


5.1 


0.5 


0.2 


Mh in MFV at large tan /? 


5 V(«o + ai)(ao + 02) 


V 50 ) 


NMFV 


62 


6.2 


2 


General 


24000 


2400 


800 



TABLE V: Summary of the 95% probability lower bound on the NP scale A (in TeV) for several 
possible flavour structures and loop suppressions. 
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